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ABSTRACT 



The Herschel Lensing Survey (HLS) takes advantage of gravitational tensing by massive galaxy clusters to sample a population of high-redshift 
galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100- 
500 ;um bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star 
formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). 
By combining our observations with LABOCA 870 yum and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 
24yum-selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with 
systematically lower infrared luminosities. This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star 
formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies 
with star formation occurring on larger physical scales. 
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1. Introduction 

Luminous (and ultra-luminous) infrared galaxies ((U)LIRGS) 
are a population of massive star-forming galaxies which con- 
tain significant amounts of dust, absorbing the ultraviolet emis- 
sion from newly formed stars within them and re-radiating this 
energy at far-infrared (FIR) wavelengths. The FIR emission of 
these sources is therefore stro ngly correlated with their level 
of star formation activity (e.g. lKennicuttlfl998l) . Spitzer MIPS 
measurements from 24-160//m bracket the peak of the dust 
emission from nearby (U)LIRGS. Therefore, MIPS observations 
have been very successful in characterizing the spectral energy 
distributions (SEDs ) and star formati on rates (SFRs) of the lo- 
cal population (e.g. lRieke et al.ll2009h . Similar measurements at 
FIR/submillimeter (submm) wavelengths are required to con- 
strain the SEDs of their higher redshift count erparts. Together 
the PACS (Poglits chetaDl20 10) and SPIRE (Griffin et al. 2010) 
instruments on-board Herschel (Pilbratt et al..20 10) provide sen- 
sitive measurements of this population at five wavelengths from 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with 
important participation from NASA. Data presented in this paper 
were analyzed using "The Herschel Interactive Processing Environment 
(HIPE)," a joint development by the Herschel Science Ground Segment 
Consortium, consisting of ESA, the NASA Herschel Science Center, 
and the HIFI, PACS, and SPIRE consortia. 



100-500 //m, sampling the rest-frame peak of the FIR emission 
out to a redshift of z = 4. 



The depth of submm maps is ultimately limited by confusion 
noise resulting from a high density of sources relative to the an- 
gular resolution of the telescope. Gravitational lensing by mas- 
sive galaxy clusters provides the only means to sample an abun- 
dant population of intrinsically faint or high-redshift infrared 
(IR) galaxies which lie below thi s limit. The Herschel Lensing 
Survey (HLS) (lEgami et al.ll2010h will target ~ 40 galaxy clus- 
ters with PACS and SPIRE to compile the first significant submm 
catalog of these galaxies. In this letter we present initial re- 
sults from our science demonstration phase observations: 8 'x 
8' PACS maps and 17 'x 17' SPIRE maps centered on the Bullet 
clus ter (1E0657-56). D etails of these observations are presented 
in Egam i et aP (12010). This is a unique target because of its 
strong lensing potential and the richness of existing ancillary 
data. We present the FIR properties of 19 MIPS 24yum-selected 
galaxies which are located behind the Bullet cluster We com- 
pare the shapes of their SEDs with those of local (U)LIRGS, and 
compare their measured IR luminosities with predictions which 
extrapolate this quantity from the observed MIPS 24 /im flux. 
Subsequent uses of MIPS in the text refer to the MIPS 24yum 
band. 
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Fig. 1. The most highly magnified galaxy in our sample, HLS 1 8, is the same LIRG first detected bv lWilson et alJ (l2008l) with AzTEC 
at 1 . 1 mm. Left panel: 1 .5'x2' thumbnails of the source at a variety of wavelengths highlight the importance of multi-wavelength 
data to guide the interpretation of the lower resolution maps. IRAC data reveal three lensed images of the galaxy indicat ed with 
white circles, although only emission from the two brightest images is detected in the Herschel maps. Right panel: Riek e et aP 
(|2009|) SED template fit to the source. The fluxes include the contribution from both of the brightest images of the galaxy, indicated 
by the two left-most circles in the thumbnails. The data points have not been de-magnified. The 24jjm flux is plotted for reference, 
although i t was not used in the fit. The best-fit modified blackbod y is shown with a dotted line and is consistent with previous 
estimates dWilson et al.ll2008l; iGonzalez et al. 1120091 iRex et al.ll2009l) . 



2. Data and source catalog 

The PACS and SPIRE observations of the Bullet cluster were 
reduced using the Herschel Interactive Processing Environment 
(HIPE). Small deviations fr om the standard pipeline are dis- 
cussed in lEgami et al.1 (l2010h . 



2.1. Source selection 

PACS and SPIRE data provide unprecedented sensitivities at 
FIR/submm wavelengths. Nonetheless, the precise identification 
of IR galaxies based solely on Herschel maps is complicated due 
to relatively large beam sizes leading to confusion noise in the 
SPIRE bands. We have therefore constructed an initial catalog 
with the positions of IR sources in the region based on secure 
(> 10-cr) detections in the higher resolution Spitzer MIPS 24yum 
map of the field. A 10-cr 24 //m threshold corresponds to the 1- 
o" error in the PACS 100/im map, our deepest Herschel map of 
the field. We select the positions for our photometry from this 
MIPS catalog in an effort to minimize the effect of Eddington 
bias which artificially boosts the flux of sources selected from 
confused submillimeter maps because of the steep underlying 
number counts. A proper statistical treatment of this effect will 
be presented in the more comprehensive analysis of this field in 
preparation. 

In order to identify a subset of these initial sources which 
lie behind the Bullet cluster, we have compiled all of the spec- 
troscopic redshift information for galaxies in the field obtained 
thus far A description of these redshift catalogs an d our method 
of a ssociation with the MIPS positions is given in iRawle et alJ 
(I2010.) . A histogram of the redshifts reveals a large group of 
Bullet cluster members at z ~ 0.3, as we ll as a smaller clust er of 
galaxies located just beyond at z ~ 0.35 ("Rawl e et al.ll201 0). We 
have therefore chosen a lower limit of z > 0.40 to select back- 
ground field galaxies. With these criteria we assemble a sub- 
catalog of 50 MlPS-selected galaxies which are spectroscopi- 
cally confirmed to be located behind the Bullet cluster. In many 
cases one SPIRE beam contains more than one MIPS galaxy. In 
these instances we are sometimes able to use the PACS resolu- 
tion to identify the dominant source of the SPIRE emission. If 
this is not possible we exclude the source from our present anal- 



ysis. Finally, we constrain our work to sources with > 3-cr detec- 
tions in at least two Herschel bands. We are thereby left with a 
sample of 15 significant background galaxies with spectroscopic 
redshifts. 

We also present the analysis of 4 additional galaxies selected 
because of photometric redshift (photo-z) estimates suggesting 
they are at Zphot > 1-5, and w hich are detected in LABOCA 
870 yum ('Johansson et al.' '2010") and AzTEC 1.1mm maps of 
the field (Wilson et al. 2008). We have calculated photo-z's for 
these sources using tw o methods: one based on IRAC colors 
(^Perez-Gonzalez et al.li 2005|) and the other on the FIR-mm col- 
ors (Hughes et al. 2002: l Aretxaga et al.ll2003h . When the two es- 
timates disagree we have chosen the value which yields the low- 
est ;i'^ fit to our SED templates. Including these galaxies in our 
analysis allows us to take advantage of the substantial correla- 
tion between SPIRE maps and longer wavelength submm/mm 
maps of the field. These results demonstrate the strength of com- 
bining such data sets to identify a potentially higher redshift 
galaxy population. The positions of our selections, along with 
corresponding MIPS and Herschel photometry can be found in 
table [T] (included in the online supplementary material). 



2.2. Photometry 

Herschel photometry for these galaxies is measured by simulta- 
neously fitting the PSF at the positions of all of the 10-cr MIPS 
sources identified in our original catalog using routines included 
in the IRAF package DAOPHOT. To ensure the maps are aligned 
to the same astrometry small offsets are applied to the PACS data 
based on the results of stacking the lOOjum map on the MIPS 
positions and to the SPIRE data based on the results of stacking 
the 250jum map on the MIPS positions. The 500 yum map is also 
corrected for contamination from th e Sunyaev-Zerdovic h effect 
based on fits to the data presented in lZemcov et alj (1201 Oh before 
performing photometry. 

In cases where multiple MIPS galaxies fall within half the 
beam size at the corresponding waveband, we approximate the 
position of the group by the average of the counterparts weighted 
by their 24jum signal-to-noise. We then take an iterative ap- 
proach; fitting at the positions of the brightest sources first. 
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Fig. 2. Left panel (a): SED template fits to HLSll, a redshift z ~ 1 LIRG. The data points have been de-magnified by a factor of 
1 .4 based on our lensing models of the foreground cluster. Although both the R09 and the CEOl templates predict the IR luminosity 
based on observed 24 //m emission to within a factor of 2, the 24//m-predicted SEDs do not provide reasonable fits to the FIR- 
mm data. Right top panels (b & c): The IR luminosity of the best-fit template as it corresponds to the local population is plotted 
as a function of the actual IR luminosity of each galaxy. The SEDs of our sample are better fit by templates corresponding to 
systematically lower luminosity local galaxies. Right bottom panels (d & e): The 24yum-predicted IR luminosity is plotted as a 
function of the actual IR luminosity of each galaxy. Both template families predict reasonable values for the luminosities of the 
galaxies in our sample, except in the cases of the ULIRGS/Hyper-LIRGS atz > 1.5. Note that HLSll &HLS18 (Fig. [B fall on the 
trend line even though their SED shapes show significant deviation in (b) & (c). 



removing them from the map, fitting at the positions of the 
next tier, and so on until we have reached the 3-cr noise level 
of the observation. This method is similar to that outlined in 
[Perez-Gonzalez et al. (2010), although we force the algorithm 
to fit the PSF at the MIPS positions rather than allowing for any 
re-centering. 



Our analysis includes LABOCA 870 //m and AzTEC 1 . 1 mm 
photometry for our sources when it is of > 3-cr significance. Four 
of our objects have LABOC A counterparts within 8" listed in the 
photometric catalog from J ohansson et al.l t2010) . In these cases 
we use the deboosted flux given in the catalog for our analysis. In 
order to obtain photometry for the remaining sources, we mea- 
sure the flux in an aperture of 40" at the MIPS position in order 
to get the total flux from the beam. Th is method is not the same 
as that used in [Johansson et al.l (l20I0i) . although it gives consis- 
tent results for the 1 7 sources presented therein. In the case of the 
four sources chosen, in part, due to bright long-wavelength emis- 
sion, we identify clear associations between the MIPS source 
and a significant, individual AzTEC source. In these cases we 
use the deboosted flux from the AzTEC catalog. For the other 
sources the AzTEC fluxes are measured at the MIPS positions 
from a PSF-convolved map. In the instances where there is a 
counterpart (within 8") in the AzTEC catalog, the difference be- 
tween the flux at the MIPS position and the deboosted point- 
source flux is < 10%. 



3. Results 

Figure [T] shows the most highly magnified galaxy in our sample, 
first r eported as a strongly lensed IRAC source (Bradac et al] 
20061: [Gonzalez et ani2009l) and independe ntly found to be the 
brightest submm/ mm source in the fi eld dWilson et al.l 120081 : 
iRex et al . 2009; Jo hansson et al.ll20IOl) . Multiple images of the 
galaxy are indicated with white circles. The right panel shows 
our SED fit to the sum of the fluxes from the brightest of the two 
images. SPIRE data indicate a 250 /im flux dens ity which is 30% 
lower than that measured in the BLAST data (IRex et al.ll2009l) . 
although the measurements are consistent within the uncertain- 
ties. The Spitzer InfraRed Spectrograph (IRS) s pectrum of this 
source confirms that it is at redshift z - 2.79 (iGonzalez et al.l 
I2OIOI) . consistent with published photo-z estimates similar to the 
ones used for the four sources in our sample which lack spec- 
troscopic redshift information. At that redshift, our model in- 
dicates a magnification of ~ 50 due to gravitational lensing by 
the foreground cluster This val ue is lower than other estimates 
presented in the literature (e.g. iGonzalez et al.l (l2OI0l) suggest 
a magnification of ~ 100). We note that this value should be 
treated as a lower limit since there is likely to be additional local 
lensing from nearby objects. We derive an observed IR lumi- 
nosity of 2.8 X 10'^ L0 based on the template SED fitting (see 
below for a detailed discussion of this SED fitting procedure). 
Adopting our estimated lower limit for the magnification of the 
source indicates an intrinsic IR luminosity of < 5 x lO" Lq 
which is c onsistent with calcii l ations previously reported in the 
literature dWilson et al.1 120081: IGonzalez et aH 120091: iRex et al.1 
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l2009t iGonzalez et alJl20Tol) . apart from the assumed magnifi- 
cation factor The detection of such an intrisically faint galaxy is 
not possible without gravitational lensing, clearly demonstrating 
the power and promise of our strategy. 

3.1. SED properties 

The galaxies in our sample span a redshift range of 0.40 < z < 
3.24. Knowing their redshifts, we can derive their apparent IR 
luminosities. We use the lensing models described in Paraficz 
et al. (in prep.) to correct for the magnification of each source 
due to gravitational lensing by the foreground cluster and cal- 
culate the galaxies' intrinsic properties. In our disc ussion of IR 
prope rties, we adopt nomenclature similar to that in lRieke et al.l 
(120091) . defining the total IR luminosity (Ltir) as the luminosity 
in the rest-frame wavelength range A - 5-1 000 /im. The red- 
shifts, magnifications, and de-magnified IR luminosities of each 
source are listed in table|2] We assume a ACDM cosmology with 
D.ni = 0.3, Q.A = 0.7, and Hq - 70km s 'Mpc ' for our calcula- 
tions. 

We use a y^ minimi zation routine to fit the SED templates 
pr esented in[R ieke et al. (|2009|) (R09), as well as those presented 
in lCharv & Elbazi (i2001 ') (CEOl). These templates are based on 
data from local galaxies and each is designated by a luminosity 
class corresponding to the total infrared luminosities of the local 
(U)LIRGS for which it provides the best fit. The left panel of 
Fig. |2] shows an example of the template fits to a galaxy in our 
sample. The solid lines show the best-fit R09 template (in blue) 
and CEOl template (in red) to the FIR-mm data, excluding the 
24/im point. The dotted lines show the respective fits to only the 
observed 24 /zm point. The SED fits to the remaining galaxies in 
our sample are shown in Figs. [3] and |4] (included in the online 
supplementary material). 

In general we find that the templates provide good fits to 
the FIR-mm data in our galaxy sample and that the IR lumi- 
nosities derived using the best-fit templates from both R09 and 
CEOl agree within a I-o" spread of 15%. We also find that the 
galaxies in our sample are best fit by templates with system- 
atically lower luminosity classes. The top right panels (b & c) 
of Fig. |2] illustrate this trend. The luminosity class of the best- 
fit template to each galaxy is plotted as a function of the ac- 
tual luminosity for the R09 templates in (b) and the CEOl tem- 
plates in (c). Although we find a large scatter in the luminosity 
classes of best-fit SEDs, both template families show the same 
general trend. This is because the SEDs of our galaxies peak at 
longer wavelengths than local galaxies of comparable luminosi- 
ties. Therefore their spectral shapes more closely resemble those 
of lower luminosity (U)LIRGS in the local universe. The result 
supports evidence that high-redshift submm galaxies are cooler 
than local galaxies with similar IR luminosities, suggesting the 
star formation within them is occurring on more extended phys- 
ical scales (e.g. Chapman et al. 2004). 

The discrepancy between the FIR SED shapes of our sources 
and those of local counterparts with similar luminosities implies 
that a blind application of the local template SED models might 
not accurately predict the Ljir and hence the SFRs of higher red- 
shift submm galaxies. However, other studies based on Herschel 
data show that the Ljir predicted by observed 24yum emission 
is in good agreement with that measur ed from the FIR data for 
galaxies at z < 1.5 (lElbaz et al.ll2010l) . The bottom right pan- 
els of Fig. |2]confirm this trend. The IR luminosity derived from 
24/im is plotted as a function of the actual IR luminosity for the 
R09 templates in (d) and the CEOl templates in (e). Although the 
R09 fits show a slightly higher deviation, both template families 



predict reasonable values for the luminosities of the galaxies in 
our sample, except in the cases of the ULIRGS/Hyper-LIRGS at 
z > 1 .5. The left panel of Fig. |2]shows a closer inspection of the 
SED fits, revealing that although the predicted IR luminosities 
are consistant with the measured values, the SED templates se- 
lected based solely on observed 24/im emission typically peak 
at shorter wavelengths, and do not provide good fits to the FIR- 
mm data. 



4. Conclusion 

We have presented the first Herschel analysis of galaxies located 
behind the Bullet cluster. We find that their colors are best fit us- 
ing templates based on local galaxies with systematically lower 
IR luminosities. This suggests that our sources are not like lo- 
cal ULIRGS in which vigorous star formation is contained in a 
compact highly dust-obscured region. Instead, they appear to be 
scaled up versions of lower luminosity local galaxies with star 
formation occurring on larger physical scales. A more compre- 
hensive analysis of the field in preparation will combine our full 
catalog of sources with LABOCA and AzTEC data to compile 
a larger sample of these galaxies. By studying their relationship 
to local star-forming galaxies we can better understand the pro- 
cesses that govern their evolution. 



Table 2. Source properties. 



ID 



Mag R09 Ltir [Lq] 



HLSOl 
HLS02 
HLS03 
HLS04 
HLS05 
HLS06 
HLS07 
HLS08 
HLS08 
HLSIO 
HLSll 
HLS12 
HLS13 
HLS14 
HLS15 
HLS16 
HLS17 
HLS18 
HLS19 



0.61 

0.72 
0.82 



4.1 + 



' 1.6;-[, 1.6;-^ 



"2.8!-], 1.8 

0.66 

0.43 

0.82 

0.60 

1.07 

3.24 

" 2.r-] 

0.57 ■ 

1.17 

0.45 

0.40 

2.79 

0.71 



3.V 



1.1 


2.6e+ll 


1.00 


4.0e+12 


1.0 


3.1e+ll 


1.0 


4.1e+ll 


1.16 


1.6e+13 


1.11 


5.0e+12 


1.1 


1.6e+ll 


1.2 


2.4e+10 


1.1 


8.3e+ll 


1.1 


3.1e+ll 


1.4 


8.3e+ll 


11.3 


3.5e+ll 


1.23 


4.2e+12 


1.1 


6.0e+10 


1.1 


4.7e+ll 


1.4 


2.5e+ll 


1.1 


2.7e+ll 


>54'' 


5.2e+ll 


1.1 


l.le+12 



■' IRAC photo-z, FIR-mm photo-z. Adopted value in bold. 

*" Lower limit due to unquantified local lensing by nearby objects. 
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Table 1. Flux densities measured from the MIPS and Herschel maps. The 500 yum map was corrected for SZ contamination before 
extracting the photometry. These measurements have not been corrected for the amplification due to gravitational lensing. 



ID 



RA 

[deg] 



Dec 

[deg] 



F1UX94 

[miy] 



Fluxioo 
[mIy] 



Fluxieo 
[mIy] 



F1UX250 

[mIy] 



F1UX350 

[mIy] 



Fluxsoo 
[mIy] 



HLSOl 


104.72828 


-55.88917 


0.40+ 0.01 






21.1+2.5 


10.0+ 3.5 


<7.8* 


HLS02 


104.46419 


-56.02209 


0.15±0.01 






10.1+2.4 


16.2+ 3.5 


14.4+ 4.0 


HLS03 


104.48783 


-56.03070 


0.39+ 0.02 






9.0+ 2.4 


9.5+ 3.4 


<7.8* 


HLS04 


104.45727 


-55.91533 


0.39± 0.01 






18.6+2.8 


13.2+ 3.6 


<7.8* 


HLS05 


104.60260 


-55.92013 


0.71± 0.02 


75.4+ 1.3 


164.4+5.1 


168.9+ 6.2 


120.0+4.1 


58.4+ 4.4 


HLS06 


104.64300 


-55.90990 


0.40± 0.01 


9.9+ 1.2 


28.4+ 2.4 


44.6+ 2.8 


39.0+ 3.5 


25.2+ 4.0 


HLS07 


104.60545 


-55.89051 


0.29± 0.01 


7.0+ 1.3 


7.8+2.1 


8.2+2.8 


<8.4* 


<7.8* 


HLS08 


104.63628 


-55.98351 


0.135± 0.003 


2.7+ 1.1 


4.6+2.1 


<6.0'' 


<8.4* 


<7.8* 


HLS09 


104.55591 


-55.87279 


0.65± 0.02 






23.2+ 2.6 


12.9+ 3.5 


<7.8* 


HLSIO 


104.54350 


-55.98005 


0.53± 0.02 


16.7+ 1.4 


22.6+2.1 


13.9+ 2.4 


<8.4* 


<7.8* 


HLSU 


104.56008 


-55.95848 


0.47± 0.01 


7.4+ 1.1 


15.1+2.1 


33.1+2.7 


36.9+ 3.5 


18.8+4.1 


HLS12 


104.62997 


-55.94386 


0.046± 0.003 


<2.4'' 


4.3+2.1 


7.8+ 2.5 


11.4+3.4 


<7.8* 


HLS13 


104.60567 


-55.94490 


0.28± 0.01 


2.9+ 1.1 


5.5+ 2.0 


10.1+2.5 


19.6+ 3.5 


<7.8* 


HLS14 


104.64689 


-55.88658 


0.136± 0.004 


3.0+ 1.2 


4.4+2.1 


< 6.0* 


< 8.4* 


<7.8* 


HLS15 


104.73730 


-55.88516 


0.16+0.02 






11.1+2.7 


15.3+3.8 


11.7+3.7 


HLS16 


104.58575 


-55.93920 


0.83± 0.02 


25.3+ 1.3 


51.3+2.5 


47.7+ 2.6 


24.0+ 3.6 


10.2+ 3.7 


HLS17 


104.64463 


-56.00850 


0.71+0.02 


24.7+ 1.2 


56.5+ 3.4 


49.2+ 3.4 


35.5+ 3.7 


17.9+ 3.7 


HLSISa" 


104.65471 


-55.95193 


0.49+ 0.01 


7.0+ 1.2 


24.5+2.1 


65.3+ 2.8 


98.6+ 3.9 


101.4+ 4.0 


HLSISb" 


104.65861 


-55.95057 


0.36+ 0.01 


3.8+ 1.1 


14.0+ 2.5 


< 6.0* 


< 8.4* 


<7.8* 


HLS19 


104.59877 


-55.87833 


1.10+0.03 






40.9+ 2.9 


21.6+4.0 


9.4+ 3.7 



" Multiple images of the highly lensed galaxy ( Gonzal ez et alj2003) . 
'' Upper limits correspond to statistical 3-cr noise levels of each map. 
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Fig. 3. SED template fits to HLS 1-HLS9. The data points have been de-magnified according to our lensing models of the foreground 
cluster. These magnification factors can be found in table |2l The solid lines show the best-fit R09 template (in blue) and CEOl 
template (in red) to the FIR-mm data, excluding the 24yum point. LABOCA 870 /^m and AzTEC 1.1 mm data are used to constrain 
these fits when the detection is of > 3-cr significance. However, only the MIPS and Herschel data are shown in the figure. The 
dotted lines show the respective fits based solely on the observed 24 //m point. Although both the R09 and the CEOl templates 
generally predict reasonable values for the IR luminosity based on observed 24jum emission, the 24/im-predicted SEDs typically 
do not provide reasonable fits to the FIR-mm data. 
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Fig. 4. SED template fits to HLS10-HLS19. The data points have been de-magnified according to our lensing models of the fore- 
ground cluster. These magnification factors can be found in table |2] The solid lines show the best-fit R09 template (in blue) and 
CEOl template (in red) to the FIR-mm data, excluding the lAjim point. LABOCA 870 jum and AzTEC 1.1 mm data are used to 
constrain these fits when the detection is of > 3-cr significance. However, only the MIPS and Herschel data are shown in the figure. 
The dotted lines show the respective fits based solely on the observed 24 //m point. Although both the R09 and the CEOl templates 
generally predict reasonable values for the IR luminosity based on observed 24 //m emission, the 24/im-predicted SEDs typically 
do not provide reasonable fits to the FIR-mm data. 



